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ABSTRACT 

Working memory is a critical building block for almost all cognitive tasks, and impairment can cause significant dis- 
ruption to daily life routines. We investigated the functional connectivity (FC) of the visuo-spatial working memory 
network in temporal lobe epilepsy and its relationship to the underlying white matter tracts emanating from the 
hippocampus. Fifty-two patients with unilateral hippocampal sclerosis (HS) (30 left) and 30 healthy controls 
underwent working memory functional MR] (fMRI) and Diffusion Tensor Imaging (DTI). Six seed regions were iden- 
tified for FC analysis; 4 within a task-positive network (left and right middle frontal gyri and superior parietal lobes), 
and 2 within a task-negative network (left and right hippocampi). FC maps were created by extracting the 
time-series of the fMRI signal in each region in each subject and were used as regressors of interest for additional 
GLM fMRI analyses. Structural connectivity (SC) corresponding to areas to which the left and right hippocampi 
were connected was determined using tractography, and a mean FA for each hippocampal SC map was calculated. 
Both left and right HS groups showed atypical FC between task-positive and task-negative networks compared to 
controls. This was characterised by co-activation of the task-positive superior parietal lobe ipsilateral to the typically 
task-negative sclerosed hippocampus. Correlational analysis revealed stronger FC between superior parietal lobe 
and ipsilateral hippocampus, was associated with worse performance in each patient group. The SC of the hippo- 
campus was associated with the intra-hemispheric FC of the superior parietal lobe, in that greater SC was associated 
with weaker parieto-frontal FC. The findings suggest that the segregation of the task-positive and task-negative FC 
networks supporting working memory in TLE is disrupted, and is associated with abnormal structural connectivity 
of the sclerosed hippocampus. Co-activation of parieto-temporal regions was associated with poorer working 
memory and this may be associated with working memory dysfunction in TLE. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 
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1. Introduction 

Whether the medial temporal lobe (MTL) is involved in working 
memoiy (WM) is still open to debate (Graham et al„ 2010; Jeneson and 
Squire, 2012) but it is a fundamental question pertaining to the conceptu- 
al understanding of how memory is organised in the brain. Working 
memoiy refers to the temporary storage and manipulation of information 
and has typically been considered dependent on frontal lobe integrity, 
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independent of the MTL Much of the evidence for this distinction relies 
on early studies which suggested that WM was unaffected in temporal 
lobe epilepsy patients (Cave and Squire, 1992; Drachman and Arbit, 
1966). However, emerging evidence suggests that specific WM processes 
are disrupted in TLE (Abrahams et al., 1999; Krauss et al., 1997a; Owen et 
al., 1996) and that the MTL is involved in both short and long term 
memory formations (Cashdollar et al., 2009; Ranganath and Blumenfeld, 
2005). Whether the disruption of WM in TLE is a result of critical MTL 
involvement in WM processes (Corcoran and Upton, 1993), or secondary 
to propagation of epileptic activity from the epileptogenic zone to 
eloquent cortex responsible for WM function (Hermann et al., 1 988) is 
still unknown. 

Previous imaging studies provide evidence for hippocampal in- 
volvement in specific WM processes, particularly in the context of 
spatial processing (Bird and Burgess, 2008) and increasing WM load 
(Axmacher et al, 2008). Our recent findings using a visuo-spatial 
WM fMRI task suggested that in TLE the ipsilateral hippocampus is 
unable to deactivate in response to increasing load, to the detriment 
of WM performance (Stretton et al., 2012). This study aimed to assess 
the functional connectivity patterns of the diseased hippocampus to 
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further understand the mechanisms of WM dysfunction in TLE. Previ- 
ous studies report functionally segregated networks subserving WM; 
task-positive networks, involving fronto-parietal activation (Owen et al., 
2005; Palacios et al., 2012) and task-negative networks, involving deac- 
tivation of the default-mode network (Palacios et al., 2012; McGill et al., 
2012), including the hippocampi (Cousijn et al., 2012a). Disruption of the 
segregation of task-positive and task-negative networks in resting-state 
fMRI has been described in Idiopathic Generalised Epilepsy (IGE) (McGill 
et al., 2012) but has not yet been assessed in a TLE cohort. 

Functional connectivity (FC) analysis of fMRI data identifies tempo- 
rally correlated regions to elucidate the network dynamics for the un- 
derlying function. Several studies have investigated resting-state FC 
networks in TLE patients (Bettus et al., 2009a; Liao et al., 2010, 2011; 
Morgan et al., 201 1 ; Waites et al., 2006) and during cognitive fMRI of 
WM (Axmacher et al., 2008; Campo et al., 2011), episodic memory, 
(Addis et al., 2007; Voets et al., 2009; Wagner et al., 2007) and language 
(Bonelli et al., 2012; Protzner and McAndrews, 2011; Vlooswijk et al., 
2010), with mixed results. There is evidence for both increased connec- 
tivity and decreased connectivity in multiple resting state networks in 
TLE compared to non-TLE partial epilepsy patients (Luo et al., 2011). 
Furthermore, with regard to lateralization of effect, decreased ipsilateral 
and increased contralateral resting connectivity of the MTL epilepto- 
genic network has been observed (Bettus et al., 2009b). In task-related 
connectivity studies, during complex scene encoding, the ipsilateral 
MTL has been shown to have reduced connectivity to the contralateral 
MTL and task-dependant eloquent cortex compared to controls (Voets 
et al., 2009). More recent work suggests that memory function in TLE 
may not only rely on hippocampal function alone, but also prefrontal 
network integrity, which shows reduced FC in cryptogenic localization- 
related epilepsy (Vlooswijk et al., 2011). 

The underlying structure of temporal and extra-temporal white 
matter has also been shown to influence default mode network con- 
nectivity in TLE (Liao et al., 201 1 ). Therefore in addition to FC, we de- 
rived a measure of structural connectivity and integrity of the white 
matter tracts surrounding the hippocampi. The primary measure of 
white matter integrity used in this study is fractional anisotropy (FA), 
which is determined by the directional magnitude of water diffusion 
in three-dimensional space. Tightly packed white matter fascicles result 
in water diffusion in a preferred direction (high FA), reflecting greater 
structural integrity. In contrast, white matter fascicles that have poor 
structural coherence will allow water to diffuse more randomly (low 
FA) reflecting white matter abnormalities. Finally, recent graph theory 
work suggests that the age of epilepsy onset is associated with lower 
temporal lobe FC, and longer duration of epilepsy is related to more 
random configuration of FC networks (van DE et al., 2009). 

To date, no studies have investigated the segregation of task-positive 
and task-negative WM networks in TLE, or their relationship to the un- 
derlying structure of the temporal lobes and clinical variables such as 
age of epilepsy onset and seizure frequency. We hypothesised: 

1) The segregation of task-positive and task-negative networks is 
disrupted in TLE. 

2) This disruption is associated with poorer WM performance. 

3) The structural connectivity of the hippocampus would influence 
FC network efficiency. 

4) An early seizure onset and frequent seizures would be associated 
with greater disruption of WM networks. 

2. Materials and methods 

2.1. Subjects 

Fifty-two patients with medically refractory TLE and unilateral 
hippocampal sclerosis (HS) (30 left HS, mean age 38 years, range 
19-54 years; 22 right HS, mean age 41 years, range 21-56 years) un- 
dergoing pre-surgical evaluation at the National Hospital for Neurology 



and Neurosurgery (NHNN) participated in this study. The study was 
approved by the NHNN and the Institute of Neurology Joint Research 
Ethics Committee, and written informed consent was obtained from 
all subjects. 

All patients underwent structural MR1 at 3 Tesla (3 T). Prolonged 
Video-EEG recordings confirmed seizure onset in the medial temporal 
lobe ipsilateral to the HS, and all patients had a normal contralateral 
hippocampus (volume in the range of 2.16-3.16 mm 3 ; T2 relaxation 
time< 118ms), based on qualitative and quantitative MR1 criteria 
(Woermann et al., 1998). All patients were taking anti-epileptic 
drugs (AEDs), were native English speakers and had undergone a 
neuropsychological evaluation as part of their presurgical investiga- 
tions. Clinical and demographic data are detailed in Table 1. We also 
recruited 30 native English speaking healthy volunteers (mean age 
37 years, range 19-64 years) without any history of neurological or 
psychiatric disease as controls. 



2.2. MR data acquisitions 

MRI studies were performed on a 3 T GE Excite HDx scanner. Stan- 
dard imaging gradients with a maximum strength of 40 mT m _1 and 
slew rate 150 T m _1 s _1 were used. Data were acquired using an 
eight-channel array head coil for reception and the body coil for 
transmission. Gradient-echo planar T2*-weighted images were ac- 
quired, providing blood oxygenation level-dependent (BOLD) con- 
trast. Each volume comprised 50 oblique axial 2.4 mm slices (with 
0.1 mm gap) covering the whole brain, with a 24-cm field of view, 
SENSE factor 2, 64x64 matrix, and an in-plane resolution of 
3.75x3.75 mm. Echo time (TE) was 25 ms, and repetition time (TR) 
was 2.5 s. Coronal Tl -weighted volumetric acquisition with 170 con- 
tiguous 1.1 mm-thick slices (256x256 matrix; in-plane resolution, 
0.9375x0.9375 mm) was used for hippocampal volume measure- 
ments. The entire length of the hippocampus was measured in the 
coronal oblique plane on alternate slices, using a manual tracing 
method (Cook et al., 1992; Lemieux et al., 2000). Volumes were calcu- 
lated by multiplying the sum of the cross-sectional area of each slice 
by the slice thickness. 



Table 1 

Clinical and demographic data. 





Controls 


Left HS 


Right HS 


N 


30 


30 


22 


Gender (male/female) 


12/18 


14/16 


5/17 


Handedness (right/left) 


27/3 


26/4 


18/4 


Age (years; mean; SD) 


37 (12.5) 


38 (9.8) 


41 (10.4) 


Right HV (cm 3 ) 


2.7 (.25) 


2.7 (.38) 


1.9 (.36)* 


Left HV (cm 3 ) 


2.6 (.31) 


1.7 (.36)* 


2.7 (.05) 


Age of onset (years) 


n/a 


14.5 (10.8) 


15.5 (11.9) 


Duration of epilepsy (years) 


n/a 


24 (14.1) 


26 (16.8) 


Complex partial seizure frequency 


n/a 


10.5 (16) 


6 (6.7) 


(per month) 








Average no. anti-epileptic drugs (range) 


n/a 


2 (1-4) 


2 (2-4) 


Anti-epileptic drugs 


n/a 


LEV = 19 


LEV=13 






CBZ = 14 


CBZ = 11 






LTG=10 


LTG = 6 






CLB = 9 


CLB = 5 






VPA = 5 


TPM = 5 






OXC = 5 


0XC = 4 






ZNS = 4 


PHT = 3 






TPM = 3 


PGB = 3 






PGB = 1 


LCS = 2 








PB = 2 








ZNS=1 








PRM = 1 



HS= Hippocampal Sclerosis; HV= Hippocampal Volume; * = Significant at p<0.05 
following one-way ANOVA; CBZ = Carbamazepine; LEV=Levetiracetam; VPA=Sodium 
Valproate; LTG = Lamotrogine; LCS = Lacosimide; TPM=Topiramate; PHT = Phenytoin; 
ZNS=Zonisimide; CLB = Clobazam; OXC=Oxcarbazepine; PGB = Pregabalin; PB = 
Phenobarbitone; PRM = Primidone. 
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2.3. Diffusion Tensor Imaging (DTI) acquisition 

In 6 of the 82 participants, (3 controls, 2 left HS and 1 right HS) DTI 
acquisition failed due to technical difficulties (loss of pulse rate monitor). 
Data were acquired using a cardiac-triggered single-shot spin-echo pla- 
nar imaging sequence with echo time = 73ms. Sets of 60 contiguous 
2.4 mm-thick axial slices were obtained covering the whole brain, with 
diffusion sensitizing gradients applied in each of 52 non-collinear di- 
rections (b value of 1200 mm 2 s _1 [5 = 21 ms],A = 29 ms using full 
gradient strength of 40 mT m _1 ) along with 6 non-diffusion weighted 
scans. Gradient directions were calculated and ordered as described 
elsewhere (Cook et al., 2007). Field of view was 24 cm, with matrix 
size of 96x96, zero filled to 128x 128 during reconstruction, giving a 
reconstructed voxel size of 1.875 mm x 1.875 mm x 2.4 mm. The first 
four scans were discarded to ensure magnetization equilibrium. 

2.4. 'Dot-Back' fMRI paradigm 

Subjects were required to monitor the locations of dots (presentation 
time: 440 ms; inter-stimulus interval: 1500 ms) within a diamond 
shaped box on the screen at a given delay from the original occurrence 
(0-, 1-, or 2-back). There were three 30-s active conditions in total 
(0-, 1-, and 2-back) presented to subjects five times in pseudorandom 
order, controlling for any order effect. In total, 15 stimuli were 
presented in each 30-s active block. Each active condition started with 
a 15-s resting baseline (the word 'Rest' appeared on the screen during 
this period). Subjects were required to move the joystick corresponding 
to the correct location of the current (0-back) dot, the location of the 
previously presented dot (1-back) or the location of the dot appearing 
two presentations before (2-back). 

2.5. Data analysis 

Imaging data were analysed with Statistical Parametric Mapping 
(SPM8) (www.fil.ion.ucl.ac.uk). The time series of each subject was 
realigned, spatially normalized to an acquisition-specific echo planar 
imaging template in Montreal Neurological Institute space, and 
smoothed with a Gaussian kernel of 8 mm full-width at half maximum. 
Trial-specific responses were modelled by convolving a delta function 
that indicated each block onset with the canonical hemodynamic 
response function to create regressors of interest. Each subject's move- 
ment parameters were included as confounds. Contrasts were defined 
to identify task-positive areas comparing high memory load ('2-back') 
against the control task ('0-back'). To identify task-negative areas, i.e. 
areas progressively deactivated with increasing task difficulty, we de- 
fined an additional contrast with values —1,-2 and —3 for the three 
conditions '0-back', '1-back' and '2-back', modelling such deactivation 
compared with the resting periods. These contrast images were subse- 
quently used to define seed regions for the FC analysis. 

2.6. Functional connectivity analysis 

Based on established working memory network nodes (Owen et 
al, 2005), four regions of interest (ROIs) were defined from the 
task-positive contrast; left and right middle frontal gyri (MFG), left 
and right superior parietal lobes (SPL). Based on our previous finding of 
disrupted hippocampal deactivation in TLE during the task (Stretton et 
al., 2012), two ROIs were defined using the task-negative contrast in 
the left and right medial temporal lobes (MTL) (Fig. 1 ). ROIs were defined 
using a combination of functional and anatomical criteria: clusters 
revealed by the combined group (left HS, right HS and controls) 
and main effect of task (thresholded at p<0.05 FWE) of the respec- 
tive regions were inclusively masked with the corresponding ana- 
tomical regions of interest (ROI) of the WFU Pick Atlas toolbox 
(Tzourio-Mazoyer et al., 2007; Maldjian et al., 2003). Regional time 
series were extracted from each ROI on the single-subject level. The 



average time series from the peak voxel within all six ROIs for each in- 
dividual were used as regressors for new general linear model fMRI 
analyses. These individual first level contrast images were then taken 
to the second-level for group analysis for each seed region. 

Analysis of variance (ANOVA) was performed with group as a factor 
to examine the main effects of each connectivity seed ROI and to high- 
light regions demonstrating more or less FC in one group compared to 
another. We report all activations at a threshold of p<0.001, uncorrected 
for multiple comparisons, if not stated otherwise. In view of our a priori 
hypothesis we performed small volume corrections (SVC) on results in- 
dicating activation within the hippocampus, defined as a 6 mm sphere 
surrounding the peak voxel in the cluster. This was adopted as we 
were testing a specific hypothesis regarding MTL activity and because 
of the low signal-to-noise ratio in the anterior temporal lobe. 

2.7. Independent component analysis 

In order to identify the task specificity of the functional connectivity 
maps, independent component analysis was carried out at a group level 
using MELODIC from FMRIB software library (http://www.fmrib.ox.ac. 
uk/fsl/) (Beckmann and Smith, 2005). A4D file of the realigned, normal- 
ized, smoothed images was created for each subject, and image data 
were prefiltered with a high-pass filter with a cut-off at 100 s. The algo- 
rithm was constrained to identify 32 independent components, com- 
mon across all subjects in each group and characterised with regard to 
their location, spatial extent and signal time course. Components were 
ranked according to their relative contribution to overall signal vari- 
ance, and no manual selection or rejection of components was carried 
out(Vollmar et al., 2011). 

2.8. Structural connectivity 

Eddy current correction of the DTI data was performed using the 
eddy correct tool in FSL (version 4.0.1) (Smith et al., 2004) and the 
Camino toolkit (version 2, release 767) was used to calculate fractional 
anisotropy (FA) images for each subject. Automated segmentation of 
the Tl -weighted volumetric scans using Freesurfer was used to identify 
the hippocampi (http://surfer.nmr.mgh.han/ard.edu/). The b = 0 image 
from the DTI scan was affinely registered to the Tl volumetric scan 
using FLIRT (http://www.fmrib.ox.ac.uk/fsl/flirt/index.html). The inverse 
transformation was applied to the segmentation to give the hippocampal 
seeds in native diffusion space. Tractography was performed from the 
hippocampal seeds employing a multi-tensor model with 1000 itera- 
tions from each voxel within the seed region as previously described 
(Winston et al., 2012). The connectivity distributions generated for 
each subject were thresholded at 0.05 and the mean FA of the areas to 
which each hippocampus was connected was determined to give a 
global measure of hippocampal structural connectivity used for subse- 
quent analyses (Fig. 1 ). Greater mean FA was interpreted as a measure 
of greater structural connectivity of the hippocampus. 

2.9. Neuropsychological assessment 

2.9.1. Digit span backwards 

The Digit span subtest from the WAIS-R (Wechsler, 1981 ) was ad- 
ministered to each participant and the digit span backwards trials 
were used as the measure of WM. The participants have to repeat 
digit strings of increasing length in the reverse order. Digit sequences 
ranged from 2 to 8 with two trials per sequence. Span size was calcu- 
lated as the highest digit sequence where both trials were successful 
(max score = 8). 

2.9.2. Gesture span 

The gesture span task (Canavan et al., 1989) requires the subject 
to copy sequences of hand gestures of increasing length up to 5 ges- 
tures. The test ends when participants make two consecutive errors 



276 



J. Stretton et al. / Neurolmage: Clinical 2 (2013) 273-281 






C) 



Fig. 1. Functional and structural connectivity seed regions and maps. A. Task-positive seed regions of the bilateral middle frontal gyri and superior parietal lobes (green, axial slice) and 
control group FC maps showing co-activation of a fronto-thalamo-parietal network between each region (yellow). B. Task-negative seed region of the medial temporal lobe (green, 
sagittal slice) and control group FC map showing co-deactivation of default mode network regions and bilateral hippocampus (blue). Images are overlaid on a single-subject Tl image 
provided by SPM8. C. Structural connectivity (SC) seed regions of the hippocampus (green, sagittal slice) and individual control SC map showing areas connected to each hippocampus 
( red). The mean fractional anisotropy (FA) value was calculated in each subject for each hippocampal SC map and used for subsequent analyses. Images are overlaid on a control subject FA 
image. 



at any given gesture set size or when the maximum of 5 gestures had 
been reached. The task was repeated with a parallel version immedi- 
ately after the first version was finished. One point was given for each 
successful trial. The mean span was calculated across trials and was 
used for the subsequent analysis (max score = 5). 

2.9.3. Motor sequences 

The motor sequence task devised by Canavan (Canavan et al., 1989) 
requires a sequence of 3 hand gestures to be repeated in the same order. 
Ten alternating sequences were administered in total. The test stopped 
after all 10 trials had been completed. One point was given for each 
successful trial. The total number of successful trials was used for the 
subsequent analysis (max score = 10). 

2.9.4. Working memory: composite score 

In order to explore the relationship between WM competence and 
neural activation patterns, and to avoid multiple comparisons, a single 
measure of WM was derived using a principal component analysis 
(PCA). Out of scanner scores for digit backwards, gesture span, and 
motor sequences as well as performance from the most demanding 
2-back fMRI condition were entered into a PCA. The derived score for 
each subject was then entered as a regressor of interest into an analysis 



of covariance (ANCOVA) in order to test for correlations between areas 
of FC and performance. Neuropsychological data were analysed using 
PAS W-vl 8 (SPSS; Chicago, IL, USA). 

3. Results 

3.1. Working memory performance 

One-way ANOVA revealed that both left HS and right HS groups 
performed significantly less well than controls in all working memory 
measures (p<0.05) apart from the right HS group in the 2 dot-back 
condition of the fMRI task (p = 0.08). Similarly, the composite working 
memory score derived from the PCA was significantly lower in patients 
than controls (p<0.001 ) (Table 2). 

3.2. Functional connectivity group maps 

Seeding from the left and right MFC and SPL, each group showed 
co-activation within a fronto-thalamo-parietal network of task-positive 
regions. Seeding from the left and right MTL, each group showed 
co-deactivation with the contralateral MTL, lateral temporal lobes and 
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Table 2 

Neuropsychological data (mean and standard deviations). 



Group (n) 


Digit Span 


Gesture 


Motor 


2 dot back % 


Principal 




backwards 


span 


sequences 




component 




C/8) 


(/6) 


C/10) 






Controls (30) 


4.7(1.1) 


3.2 (.6) 


6.4 (1.9) 


70 (20) 


.67 (.86) 


Left HS (30) 


3.6(1.1)* 


2.5 (.7)* 


4.2 (1.9)* 


51 (22.9)* 


-.48 (.84)* 


Right HS (22) 


3.3 (.9)* 


2.6 (.8)* 


4.1 (2.3)* 


56 (22.9) 


-.40 (.91)* 



HS = Hippocampal Sclerosis; * = One-way ANOVA results, significant p< 0.05 compared to 
controls. 



default mode regions including the precuneus, anterior cingulate cortex 
(ACC) and medial frontal gyrus (Fig. 1). 



3.3. Functional connectivity group differences 

There was no significant reduction in connectivity for any of the 
six seed regions in the left and right HS groups compared to controls 
but increased connectivity was observed (see Table 3). 



3.3.1. Middle frontal gyri 

The left HS group showed increased co-activation of the left MFG, 
with the left hippocampus (-20, -10, -14, p<0.05 FWE SVC) 
(Fig. 2). The right HS group showed increased co-activation of the left 
MFC with the right ACC (Cashdollar et al., 2009; Bonelli et al, 2012; 
Palacios et al., 2012). 



3.3.2. Superior parietal lobes 

The left HS group showed increased co-activation of the left and right 
SPL, with the left hippocampus (-36,-18, -16, p<0.05 FWE SVC). In 
the right HS group, the right SPL had increased co-activation with the 
right hippocampus (36, - 12, -20, p<0.05 SVC) (Fig. 2). 



3.3.3. Medial temporal lobes 

For both the left and right MTL seed regions, both left and right HS 
groups showed increased contralateral MTL co-deactivation compared 
to controls (p<0.05 FWE SVC). 



3.4. Independent component analysis 

The 32 independent components identified by independent compo- 
nent analysis described 87% of the total functional MRI signal variance 
in controls, 80% in left HS and 82% in right HS. In controls, the first com- 
ponent (explaining 6% of the total signal variance) was located in the 
bilateral MFG, left central region, bilateral SPL and anti-correlation of 
the precuneus, representing the working memory network, motor 
response and the default mode network respectively.. This component's 
signal changes were time-locked to the task timing (frequency 1/45 s), 
and the response showed a strong correlation with the cognitive load 
from the task, with higher amplitude in more difficult conditions 
(Fig. 3A). A similar component was identified in the right HS group 
(explaining 7.7% of total variance) but was less well modulated by 
task difficulty (Fig. 3B). In the left HS group, a disproportionately large 
component (explaining 35.6% of the total variance) compromising the 
bilateral MFG, left central region and to a lesser extent bilateral SPL rep- 
resented a less extensive working memory network with no evidence of 
anti-correlation with the default mode network (Fig. 3C). This could 
represent a ceiling effect in the left HS group reflected by the relatively 
equal amplitude across all three conditions. 

3.5. Regression analyses 
3.5.1. WM performance 

3.5.1.1. Middle frontal gyri. In controls, greater WM performance was 
associated with co-activation of the contralateral MFG for both left 
and right seeds. Poorer performance was associated with co-activation 
of the precuneus for both left and right MFG seeds (Fig. 4). There was 
no correlation with performance in either the left or right HS groups 
for either MFG seed region. 

3.5.1.2. Superior parietal lobes. Co-activation of the left and right SPL 
with the precuneus was associated with poorer performance in the 
control group. In the left HS group, poorer performance was associated 
with increased co-activation of the left hippocampus with the left SPL, 
and of both hippocampi with the right SPL. In the right HS group, poorer 
performance was associated with co-activation of the right SPL and 
right parahippocampal gyrus (Fig. 4). 



Table 3 

fMRI peak activation group differences for each FC seed region. 



Seed region 


Group/interaction 


Z-score 


p-value 


Peak coordinates (x, y, z) 
in MN1 space 


Anatomical region 


Left MFG 


Left HS> Controls 


3.89 


0.001 


-56,22, 20 


L. inferior frontal gyrus 






3.83 


0.001 


-10, -2, 50 


L. suppl. motor area 






3.24 


0.05* 


-20, -10, -14 


L. hippocampus 




Right HS> Controls 


3.79 


0.001 


8, 28, 16 


R. ant. cingulate cortex 


Right MFG 


Left HS> Controls 












Right HS> Controls 










Left SPL 


Left HS> Controls 


3.81 


0.001 


-40, -22, 20 


L. insula 






3.59 


0.001 


-22, -2, -20 


L. amygdala 
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Fig. 2. Functional connectivity group differences. Compared to controls, left HS patients showed greater co-activation of the left superior parietal lobe (SPL) and left middle frontal gyrus 
(MFG) with the left hippocampus, while right HS patients showed greater co-activation of the right SPL and right hippocampus. HS = Hippocampal Sclerosis; L=Left; R= Right. Images 
are overlaid on a single-subject Tl image provided by SPM8. 



3.5. 3 .3. Medial temporal lobes. There was no significant correlation with 
performance in any group with either the left or right MIL seed regions. 

3.5.2. Structural connectivity 

One-way ANOVA revealed left HS patients to have significantly 
lower FA values in both the left (p<0.0001) and right (p<0.01) hippo- 
campal SC maps compared to controls. Right HS patients showed re- 
duced FA only in the right hippocampus SC map (p<0.005). 

Inline Supplementary Table SI can be found online at http:// 
dx.doi.org/1 0.101 6/j.nicl.201 3.01 .009. 

In all three groups, the FA of the left and right hippocampus SC 
maps correlated with FC from the SPL seeds, but not with the MFC 
seeds. 

In controls, lower right hippocampal SC was associated with 
co-activation of the right SPL and right inferior frontal gyrus. Lower 
left hippocampal SC was associated with greater co-activation of the 
left SPL and bilateral MFG (p<0.01 uncorrected). 

In the left HS group, lower left hippocampal SC was associated with 
co-activation of the left SPL and the (task-negative) left medial frontal 
gyrus (p<0.01 uncorrected). However, unlike controls, greater SC of 
the right hippocampus was associated with greater co-activation of 
the right SPL and right MFG. 

In the right HS group, lower right hippocampal SC was associated 
with co-activation of the right SPL and the (task-negative) bilateral 



medial frontal gyrus. Similar to controls, lower left hippocampal SC 
was associated with co-activation of the left SPL and the left superior 
frontal gyrus. 



3.6. Correlation with clinical variables 

In the left HS group, younger age of onset was associated with in- 
creased co-activation of the task-negative left MTL and task-positive 
right MFG. In addition, greater seizure frequency was associated 
with co-activation of the task-negative left MTL and task-positive 
left inferior frontal gyrus in the left HS group. There was no significant 
correlation with hippocampal volume in any group for any seed 
region. 



4. Discussion 

4.1. Summary of main findings 

In TLE there was a disrupted segregation of task-positive and 
task-negative networks, and this was associated with poorer WM. In 
addition, the SC of the hippocampus was associated with the FC of 
the parietal seed regions. 
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Fig. 3. Independent component analysis results. The first component of the controls (A) shows task-modulated fronto-parietal network connections with anti-correlation of the 
precuneus. The right HS group (B) follow a similar pattern but with less modulatory effect of the task. The left HS group (C) first component compromised the bilateral MFG, 
left central region and to a lesser extent bilateral SPL modulated equally by each condition. HS = Hippocampal Sclerosis. 
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Fig. 4. Functional connectivity; correlation with working memory performance. Co-activation of the left SPL and left hippocampus, and right SPL with the right parahippocampal 
gyrus, was associated with poor performance in the left and right HS groups respectively. In controls, there was co-activation of the task-positive seeds with the task-negative 
precuneus. HS = Hippocampal Sclerosis; L = Left; R= Right. Images are overlaid on a single-subject Tl image. 



4.2. Comparison with previous literature 

Both left and right HS patients performed less well than controls 
across all measures of WM. This is in keeping with other evidence 
suggesting WM impairments in TLE (Abrahams et al., 1999; Owen 
et al., 2005; Axmacher et al., 2007; Black et al., 2010; Grippo et al., 
1996; Krauss et al., 1997b; Wagner et al., 2009) implying a critical 
role for the hippocampus in WM functions. Previous connectivity 
studies of WM in TLE reinforce these findings. In a recent study, 
Campo and colleagues (Campo et al., 201 1 ) compared dynamic causal 
models extracted from MEG recordings during verbal WM encoding 
in 1 1 left HS patients and 1 1 healthy volunteers. The controls performed 
significantly better than the left HS group and model comparison re- 
vealed a bilateral bi-directional model including frontal and temporal 
nodes yielded the most convincing representation of verbal WM. 
Comparing this model between groups, the left HS patients showed a 
reduced ipsilateral backward connection from the left hippocampus to 
the left inferior temporal cortex. Contralaterally, patients showed sig- 
nificantly increased forward connections from the right hippocampus 
to the right inferior frontal gyrus (IFG) compared to controls, and back- 
ward connections from the right IFG and right hippocampus which 
were associated with poorer performance (Campo et al., 2011). Our 
results add to this evidence, suggesting that the FC networks 
subserving WM in TLE are disrupted both local to and remote from 
the epileptogenic zone. 

We identified two networks associated with the WM fMRI task. 
One network comprised the task-positive regions typically associated 
with WM, namely the MFC and SPL. The fronto-parietal task-positive 
network relates to the short-term storage and manipulation of infor- 
mation needed for goal-orientated behaviour that underpins WM 
processes. The other task-negative network comprised areas of the 



default mode network (DMN) that deactivate during cognitive effort 
(Raichle et al., 2001), including the hippocampi. Deactivation of the 
hippocampus during WM has previously been observed in schizo- 
phrenia (Meyer-Lindenberg et al., 2001) and healthy volunteers 
(Astur and Constable, 2004; Astur et al., 2005; Cousijn et al., 2012b; 
Hampson et al., 2006). In controls, these networks were functionally 
segregated, in line with previous evidence (Fair et al., 2007). Both 
left and right HS groups, however, showed disrupted segregation of 
the task-positive and task-negative networks, with greater ipsilateral 
co-activation of the SPL and MTL relative to controls. Disrupted segre- 
gation of task-positive and task-negative networks has been reported 
in traumatic brain injury (Palacios et al., 2012), autism (Rudie et al., 
2012) and other epilepsy syndromes (McGill et al., 2012; Vollmar et 
al., 2011). 

Using the same fMRI paradigm, patients with Juvenile Myoclonic Ep- 
ilepsy (JME) have shown increased co-activation of the motor cortex 
and the fronto-parietal network as task demands increase. Moreover, 
the default mode regions remained relatively more active during the 
task compared to controls (Vollmar et al., 201 1 ). In ICE, using resting- 
state fMRI, default mode network abnormalities were observed showing 
co-activation between positive-frontal and negative-parietal regions 
compared to controls (McGill et al., 2012). 

Our study underlines the importance of network segregation 
utilising task-related FC. WM performance negatively correlated 
with co-activation of the task-positive seeds and the task-negative re- 
gion of the precuneus in controls. In TLE, co-activation of ipsilateral 
parieto-temporal regions was associated with poorer WM, suggesting 
that dysfunction of the sclerotic hippocampus may impair the segrega- 
tion of WM networks in these patients. The group ICA revealed 
task-specific WM network components in our controls and right HS 
group and to a lesser extent in the left HS group. In the left HS group, 
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the first component explained almost 36% of the total signal variance, 
reflecting the heterogeneity of connectivity within the left HS group. 
While this could influence the interpretation of the FC results, perhaps 
reflecting reduced task compliance, it is noted that there was no signif- 
icant difference in task performance between the left and right HS 
groups. Furthermore, the ipsilateral parieto-temporal increases and re- 
lated correlations with poor performance are also seen in the right HS 
group. However, it is of note that there was no correlation between me- 
dial temporal lobe FC and performance in any group, only from the 
task-positive superior parietal lobe. Further measures of effective con- 
nectivity to infer causality and direction of effect would be needed to 
explore this result. Finally, in the left HS group younger age of onset 
and more frequent seizures were associated with more marked disrup- 
tion of networks. This is in keeping with findings that earlier insult and 
higher seizure frequency are more disruptive to functional network ar- 
chitecture (Janszky et al., 2003; Kaaden etal., 2011). 

Thus, evidence from both the generalised and focal (TLE) epilepsy 
studies indicate that disrupted segregation of functional networks 
could be a common phenomenon of epilepsy, with syndrome-specific 
effects depending on the epileptogenic networks involved in each 
syndrome. 

We also report a novel observation that the SC of the hippocampus, 
independent of task, is associated with the segregation of FC networks. 
The SC of the hippocampus was not associated with connectivity of the 
frontal lobe seeds and only showed a relationship with the parietal seed 
regions. Lower FA values of the left and right hippocampus SC maps 
were associated with greater ipsilateral parieto-frontal connectivity in 
our control group. Interestingly, in the left HS group we observed an op- 
posite effect in the contralateral hippocampus, where greater FA of the 
right hippocampus SC map was associated with increased right-sided 
parieto-frontal connectivity. This could be interpreted as structural 
compensation of the contralateral hippocampus which has previously 
been reported at a functional level during verbal WM (Campo et al., 
201 1 ). In keeping with previous evidence, the FA of the ipsilateral hip- 
pocampus SC map in both patient groups was significantly lower than 
controls (Focke et al., 2008). However in patients lower FA of the ipsilat- 
eral hippocampus SC map was associated with co-activation of the ipsilat- 
eral task-positive SPL with task-negative DMN, again implying weaker 
segregation of WM networks with more severe pathology. 

White matter integrity has previously been shown to influence 
resting-state FC in TLE (Liao et al., 201 0) and extra-temporal regions, in- 
cluding the parietal lobe have been shown to be atrophied in TLE (Riley 
et al., 2010; Keller and Roberts, 2008), suggesting that white matter 
degradation associated with chronic uncontrolled seizure can disrupt 
the functional connectivity of regions local to and remote from the 
epileptogenic temporal lobe. These results taken together with our 
findings suggest that the temporo-parietal connectivity network 
may be as much at risk in drug-resistant TLE as frontal lobe networks. 

4.3. Strengths and limitations 

A major strength of our study is the large homogenous group of HS 
patients. Using patients with unilateral pathology in a bilateral task 
allowed us to identify laterality effects regarding the ipsilateral archi- 
tecture of WM networks in relation to hippocampal pathology. 
Assessing functional connectivity derived from an average time series 
across a cognitive task as opposed to resting state data can limit 
the interpretation of the results. As our block design includes three 
'active' conditions and resting periods the observed correlations may 
be driven, for example, by the rest condition and not by the task itself. 
However, the ICA identified task modulated functional connectivity 
networks involving the task-positive nodes we used. To understand 
the direction of effects, methods assessing effective connectivity such 
as dynamic causal modelling or psycho-physiological interactions 
could infer the causal influence between seeds that would facilitate 
our understanding of the neural networks supporting WM. Both patient 



groups performed significantly less well on the motor sequence task 
and gesture span task than the controls, which could have impacted 
on the ability to manoeuvre the in-scanner joystick. However the joy- 
stick response inside the scanner is equivalent to a gesture span of 
one, and this level of ability was reached by each participant in this 
study therefore it is unlikely that these results impacted on the use of 
the joystick. It is also important to note that although we have found 
similar results to those in other epilepsy syndromes regarding the seg- 
regation of positive and negative networks, there are distinct methodo- 
logical differences between these studies regarding the use of resting- 
state versus task-based connectivity measures. The impact of AEDs on 
FC is yet to be determined, and may significantly confound our results. 
The effect of Topiramate on cognition is well documented (Szaflarski 
and Allendorfer, 2012) and it has recently been shown that lamotrigine 
can inhibit the effective connectivity of motor circuits and increase the 
connectivity of prefrontal circuitry (Li et al., 201 1 ). 

5. Conclusion 

We identified disrupted segregation of the task-positive and task- 
negative networks supporting the critical cognitive function of working 
memory in TLE with HS. Co-activation of ipsilateral parieto-temporal 
regions was associated with poorer performance in both left and right 
HS patients and this may be associated with WM dysfunction in TLE. 
This study supports the notion of hippocampal dependent working 
memory and provides further evidence of remote effects associated 
with drug-resistant TLE from both functional and structural perspec- 
tives. Further work is needed to identify the causal influence of this 
disruption and whether this effect is specific to working memory or re- 
producible with other cognitive fMRI paradigms. Finally it will be essen- 
tial to evaluate the impact of temporal lobe surgery on the functional 
segregation of cognitive networks. 
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